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ABSTRACT 

. We use weak gravitational lensing to analyse the dark matter halos around satellite 

H ■ galaxies in galaxy groups in the CFHTLenS dataset. This dataset is derived from the 

. . .' CFHTLS-Wide survey, and encompasses 154 deg^ of high-quality shape data. Using 

the photometric redshifts, we divide the sample of lens galaxies with stellar masses 
in the range IQ^Mq to into those likely to lie in high-density environments 

(HDE) and those likely to lie in low-density environments (LDE). Through compari- 
son with galaxy catalogues extracted from the Millennium Simulation, we show that 
the sample of HDE galaxies should primarily (~ 61%) consist of satellite galaxies 
in groups, while the sample of LDE galaxies should consist of mostly (~ 87%) non- 
satellite (field and central) galaxies. Comparing the lensing signals around samples 
of HDE and LDE galaxies matched in stellar mass, the lensing signal around HDE 
galaxies clearly shows a positive contribution from their host groups on their lensing 
signals at radii of ^ 500-1000 kpc, the typical separation between satellites and group 
centres. More importantly, the subhalos of HDE galaxies are less massive than those 
around LDE galaxies by a factor 0.65 ± 0.12, significant at the 2.9a level. A natural 
explanation is that the halos of satellite galaxies are stripped through tidal effects 
in the group environment. Our results are consistent with a typical tidal truncation 
radius of ^ 40 kpc. 
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1 INTRODUCTION 

In the standard picture of hierarchical structure formation, 
larger dark matter halos are built up through the accre- 
tion, stripping, and mergers of smaller halos. At the ex- 
tremes of the halo mass spectrum, namely isolated field 
galaxies and galaxy clusters, we have a relatively good pic- 
ture of how the mass within these structures is organized. 
For isolated galaxies, most of the mass is contained within a 
halo of dark matter, as confirmed b y galaxy-galaxy weak 



2OOII: iParker et all l2005l: iMandelbaum et all I2OO6I: 



gravitational lensin g measurements i Brainerd et al. 19961: 



Hudson et ai] Il99^: iGuzik fc Seha^ l2002l: iHoekstra et all 



20031. l2004':'Mandelba um et al.ll200alvan Uitert et al.ll2oTa 

Vclandcr ct al. 2012). Simulations have shown that the 
shape of thi s halo can be well-app roximated by an NFW den- 
sity profile (iNavarro et al. 1997), wh i ch ha s been confirmed 
observationally llKleinheinrich et al.l l2003l : iHoekstra et al.l 
I2OO4! : 'Mande lbaum et al.ll200 8h. In galaxy clusters, most of 
the mass also seems to lie within an NFW dark matter halo, 
with the constituent gala^xies con tributing only small per- 
turba tions to the density profile (|Mandelbaum et EtUbOOSl . 
I2OO8I ). Gravitational lensing measurements have shown that 
the halos around individual galaxies within clusters are 
significantly smaller than the halos around comparably- 
luminous field galaxies, and this effect is mo re extreme with 
galaxies closer to the centres of clusters l|Limousin et al.l 
I2OO7I : iNataraian et al.|[2009,l . 

However, between the extremes of field galaxies and rich 
clusters, the picture is less clear. Since multiple galaxies 
must merge together to eventually form clusters, at some 
point the mass in the galaxies' individual halos must mi- 
grate into a shared halo. This process most likely occurs 
through tidal stripping: when two galaxies pass near each 
other, the particles in the halo of the less massive galaxy 
will tend to be "stripped" from it and thus join the more 
massive galaxy's halo. This effect has been demonstrated 
in va r ious N-body dark matter simulations IHay a shi et al., 
12004 iKazantzidis et all 12004 ISpringel et ail |2008| ). Tidal 
stripping is also expected to remove hot gas from less 
massive galaxies, which will have the effect of cutting off 
their supply of cold gas and quenching t heir star formation 
in a process known as "stra ngulation" l|Tinslev fc LarsonI 
1 19791 : iBalogh fc Morris! I2OO0I ). Galaxies in dense environ- 
ments are known to be significantly redder on average 
than field galaxies IIDres slcr 1980; Butcher fc Oemler 1984; 
iMoore et all 19961 : iBalogh et al.lll999l . l2004i '). and tidal strip- 
ping may contribute to the quenching of star formation, so 
there is a strong motivation to understand the mechanics 
and timing of tidal stripp ing (|van den Bosch et al.l I2OO8I : 
iKawata fc MulchaevI |2008| ). It remains unclear, however, 
whether this process is rapid or gradual. This question can 
in part be investigated through analysis of the group and 
cluster scales on which tidal stripping can be observed to 
occur. 

In this paper, we focus on galaxy groups, an intermedi- 
ate mass scale between field galaxies and clusters (typically 
structures in the mass range 10^^ Mq < Mhaio < 10^* Mq 
are considered groups, and more massive structures are 
considered clusters). Weak gravitational lensing provides 
the only practical tool to measure the density profiles and 
masses of dark matter halos around satellite galaxies within 
groups. Lensing analyses of groups (jHoekstra et al.l 



Johnston et al. 20071: Hamana et al. 20091 : Leauthaud et al.l 
20ld : lFord et al.ll2012l ~ have shown that the group lensing 
signal can be measured and is on average consistent with 
an NFW density profile. However, it is unclear how much 
of this signal results from a central ha lo, and how muc h 
is due to the contributions of satellites (jGillis et al.ll20l3 V 
As such, it is necessary to measure the lensing signals 
around satellites themselves to get a full picture of the 
mass distribution. Only limited wo rk has been done in the 
group regime to date. For example, ISuvu fc Halkolal (|2010l ) 
studied a strong-lensing system and determined that tidal 
stripping did seem to occur around the satellite studied, 
which lies in a group of mass on the order of 10^^ M©. 
While this result is promising, a broader base of data will 
be needed to develop a general understanding of the dark 
matter properties of satellite galaxies in galaxy groups. 

In principle, it is possible to study tidal strip- 

! iing using spectroscopically-derived gr oup catalogues 
Pastor Mira et al.l [ioTlf : iGilhs et"aLl I2OI2I '). Spectroscopic 
data allow one to more accurately identify whether a given 
galaxy is a "central" or a "satellite'. By using velocity infor- 
mation in addition to projected separation, it is also possi- 
ble to assess statistically whether a satellite is falling in for 
the first time, or has passe d pericentre and hence is tidally 
stripped (|Oman et al.ll20l"3 ). However, such analyses require 
a large galaxy sample with both spectra and deep imag- 
ing data, and such data are expensive to acquire. However, 
photometric redshifts are often available alongside imaging 
data. Due to their large uncertainties, photometric redshifts 
have the drawback that groups are difficult to detect, and 
group-central galaxies are very difficult to identify. This then 
calls for a statistical approach, calibrated by simulations, to 
simultaneously fit both the satellites' and groups' contribu- 
tions to the stacked lensing signal. This is the approach we 

will take in this paper. 

A similar approach was taken by IMandelbaum et al.l 

l|2006l ), who investigated a large selection of galaxies from 
the SDSS, selected by environment. They tested for the 
presence of stripping by fitting models for the lensing sig- 
nal to the galaxies in high-density environments (HDE) 
and low-density environments (LDE). The authors found 
no significant evidence of tidal stripping, but did not com- 
ple tely rule it out, either. This was also later attempted 
bv Ivan Uitert et al.l (|201ll ). who used the overlap between 
RCS2 and SDSS, but they were similarly unable to get a 
clear detection of tidal stripping. Here we will use data from 
the CFHTLenS collaboration, which is significantly deeper 
than the SDSS, and hence should provide a stronger lensing 
signal. We also apply a new environment estimator, which 
is tuned to work for photometric redshifts (see Section 
and a modified halo model designed to work with this envi- 
ronment estimator (see Section [3|. 

In Section [5] of this paper, we discuss the datasets used 
in the analysis and the algorithm for estimating galaxy envi- 
ronments. In Section (3] we detail the models for the lensing 
signals and the procedure used to fit the models to the mea- 
sured signals. In Section (4] we present the results of the 
analysis and discuss possible sources of error. We conclude 
in Section O 

For consistency with the Millenium Simulation, we use 
the following cosmological parameters: = 73 km s~^ 
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Mpc"\ flm = 0.25, = 0.75, and fib = 0.045. All stated 
magnitudes are in the AB system. Since there is no clear di- 
vision between galaxy groups and galaxy clusters, we use the 
terminology "galaxy groups" throughout this paper, even 
though some of the structures we refer to as such would be 
more commonly deemed clusters. When masses are quoted 
in this paper, M is used to refer to the total (halo + stellar) 
mass of a galaxy or group, and m is used to refer to the 
stellar mass of a galaxy, unless otherwise specified. When 
radial measurements are used in this paper, R refers to a 
projected, 2D proper distance. All masses are in units of 
Mq unless otherwise specified. 



rough agreement wi th deeper data suc h as WIRDS, which 
includes NIR filters l|Bielbv et al.ll2012l '). up to z = 0.8. 

Since we perform a differential measurement between 
samples, an overall bias in the stellar masses would not af- 
fect our results. It is possible, however, for a relative bias 
in the stellar mass estimates of red and blue galaxies to im- 
pact our results. This possibility is investigated further in 
Section liXn 

For this paper, we use take all unmasked galaxies with 
photometric redshifts in the range 0.2 < 2:phot < 0.8 as lens 
candidates. We divide these into HDE and LDE samples as 
described next. 



2 DATA AND SIMULATIONS 

In this section, we discuss the datasets used and how galaxies 
are selected for the HDE and LDE samples. In Section [2. 11 
we discuss the CFHTLenS survey, from which we draw our 
data. In Section [2.21 we discuss the algorithm used to esti- 
mate the local density around galaxies in the sample. Sec- 
tion 12.31 describes how the sample is divided into matched 
high-density and low-density subsamples, and presents the 
statistics of the galaxies in the HDE and LDE samples. In 
Section [2.41 we discuss the simulations we have run to test 
our methods and the statistics of the galaxy samples within 
these simulations. 



2.1 Observations 

CFH TLenS is a 1 54 d eg^ survey (125 deg^ after mask- 
ing) jErben et all |2012| ). based on the Wide component 
of the Canada - Franc e-Hawaii Telescope Legacy Survey 
mans et al.ll2012l '). which was observed in the period 
from March 22nd, 2003 to No vember 1st, 2008, using the 
MegaCam instrument tBoulade et al.l [20031 ). It consists of 
deep, sub-arcsecond, optical data in the u*g'r'i'z' filters. 
CFHTLS-Wide observations were carried out in four high- 
galactic-latitude patches: 

• Wl: 72 pointings; RA=02'n8"'00', Dec=- 07^*00"" 00' 

• W2: 25 pointings; RA=08''54'"00', Dec=-04'*15"'00' 

• W3: 49 pointings; RA=14''17'"54=, Dec=+54'^30""31' 

• W4: 25 pointings; RA=22''13'"18', Dec=-h01'*19"'00'. 

Shapes have been measured with the LensFit shape 
measurement algorithm for galaxies with i' < 24.7 
(|Miller et al.i r2012'). giving an effective galaxy density of 11 
sources/arcmin^ in th e redshift range 0.2 < Zphot < 1-3 
|h^ mans et al.ll2012l ). Photometric redshifts are available 
for the entire s urvey, with a typical red shift uncertainty of 
~ 0.04(1 + z) iHildebrandt et al.1 120121 '). We use all fields 
in the survey, not simply those that passed the systematic s 
tests for cosmic shear measurements (|Hevmans et al.ll2012l ). 
It has been demonstrated that fields with systematics that 
may affect cosmic shear have no effect g alaxy-galaxy lens- 
ing measurements (|Velander et al.ll2012l ). and the analysis 
in this paper requires as many lens-source pairs as possible. 

We use the st ellar mass estimates described by 
IVelander et al.l (|2012l '). obtained by fitting spectral energy 
distribution (SED ) templates, following the method of 
lllbert et al.1 (|201(]| '). These stellar masses were found to be in 



2.2 Determining Environment: The P3 Algorithm 

It is not a trivial matter to determine which galaxies are 
members of groups. Even when spectroscopic redshifts are 
available, the peculiar velocities of galaxies make it impos- 
sible to de termine the membershi ps of groups with absolute 
certainty jRobotham et al^l2011^ . When only photometric 
redshifts are available, the best we can do is to select galax- 
ies that are likely to be members of groups. To do so, we use 
a modified version of the Photo- z Probability Peaks (P3) 
algorithm l|Gillis fc HudsonlbOll] '). The P3 algorithm gen- 
erates a 3-D density field by smoothing the distribution of 
galaxies in the redshift direction according to the probabil- 
ity distribution function of their photometric redshifts. The 
algorithm identifies peaks in the pseudo-three-dimensional 
field with group centres. Here we do not use the group cen- 
tres, but rather use the entire P3 density field to identify 
overdense regions. Rather than use the local P3 overdensity 
itself, we restrict ourselves to regions in which we have high 
confidence in the overdensity, and instead use the signal-to- 
noise (S/N) of the local overdensity, under the assumption 
that galaxies in overdense regions of space are more likely 
to be in groups than galaxies in underdense regions. 

We now briefly review the technical details of the P3 
algorithm. To determine the S/N of a given test galaxy, 
the P3 algorithm compares the density of gala:xies within 
a circular aperture {R = 0.5 Mpc) surrounding each test 
galaxy to the density of galaxies within a larger annulus 
(i?inncr = 1 Mpc, -Router = 3 Mpc) surrouuding each test 
galaxy (to approximate the background density). The con- 
tribution of each galaxy to this measurement is weighted by 
the probability that this galaxy lies at the same redshift as 
the test galaxy (by taking the integral of the photo-z prob- 
ability distribution function over a thin redshift slice). This 
gives the overdensity: 



5 = 



pap pannu 



(1) 



where pap and pannu are the weighted densities of galaxies 
within the aperture and annulus surrounding the test galaxy, 
respectively. This value can take the range — 1 < 5 < oo, 
which negative values corresponding to regions less dense 
than the background density, and positive values to over- 
dense regions. We then estimate the noise in this value by 
assuming a Poisson distribution for galaxies; 



f^Poisson 




+ 



/ Pannu \ 
\ T^annu / 



(2) 
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where riap and riannu are the numbers of galaxies in the aper- 
ture and annulus respectively with more than a threshold 
weighlQ. From this, we calculate the S/N = (5/(jpoisson for 
each test galaxy. Note that this S/N can take negative val- 
ues, assuming 5 is negative. The distribution of galaxies' 
S/Ns that results from this calculation depends on the choice 
of threshold weight used, so our choices of S/N limits are not 
universally applicable. We picked limits of S/N > 2 for the 
high-density sample and S /N < for the low-density sample 
based on an analysis of the simulated galaxy catalogues to 
maximize the expected signal for tidal stripping0 

Since this environment estimator provides us with 
galaxy samples biased to lie in high- and low-density en- 
vironments, we cannot use the standard halo model (eg. 
iMandelbaum et al.1 l2006l : IVelander etHI |2012| ) for fitting 
our lensing signals. Instead, the models we use are calibrated 
from simulations and are detailed in Section [3] 

2.3 Galaxy Matching 

2.3.1 Matching Algorithm 

We use the S/N values obtained for each galaxies in Sec- 
tion l2.2l to form two samples of galaxies from the catalogues. 
As we cannot ensure that a pair of random samples of galax- 
ies in high- and low-density environments will have the same 
distribution of stellar mass and redshift as each other, we 
perform a matching between galaxies with S/N > 2 and 
galaxies with S/N < as follows: 

(i) For each galaxy with S/N > 2, we search through all 
galaxies with S/N < within the same pointing. 

(ii) For each S/N < galaxy, if its stellar mass differs 
from the stellar mass of the S/N > galaxy by more than 
20% of the latter's mass, we exclude it as a possible match. 

(iii) For each remaining S/N < galaxy, we calculate a 
quality-of-match value: 



+(10(logmH-logmL))' (3) 

where zh and zl are the redshifts of the S/N > 2 and S/N 
< galaxies, respectively, and mn and ttil are their stellar 
masses. This form significantly prioritizes a match in mass 
over redshift. 

(iv) We select the four S/N < galaxies with the lowest d 
values as matches for this S /N > 2 galaxy. If there are fewer 
than four match candidates, we assign them all as matches. 

(v) Assuming at least one match was found for it, we add 
this S/N > 2 galaxy to the HDE sample, and we set its 

^ We use a threshold weight here of a > 0.001% chance of lying 
within a redshift of 0.01 of the test galaxy. 

^ In a rough approximation, the expected signal-to-noise 
of a stripping measurement is proportional to (/sat, HDE — 

/sat,LDE)\/ ^HDE + ^LDE' "here /sat, HDE and /sat,LDE are the 
fractions of satellites in the HDE and LDE samples respectively, 
and A^HDE and A'lde are the number counts of galaxies in the 
HDE and LDE samples respectively. We calculated this value for 
various S/N cuts, and the combination of S/N > 2 for the HDE 
sample and S/N < for the LDE sample provided the best ex- 
pected signal-to-noise for a stripping measurement. 
^ Matching only within the same pointing is done to conserve 
computational time. 
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Table 1. Statistics of galaxies in various stellar mass bins in the 
CFHTLenS survey, as a function of environment, z is the mean 
redshift of the bin. /^gj is the fraction of galaxies that are red, 
and /blue is the fraction that are blue, determined by the best- 
fit photo metric template s and defined in the same manner as 
by .Velander et al.l l|2012h . Fractions do not add to unity as not 
all ga laxies are classified as "red" or "blue." See I Velander et al.l 
l|2012l ) for further explanation. All average values and fractions 
assume galaxies are weighted by their stellar masses. 

weight equal the number of matches we found. (This weight 
is later applied when we stack lensing signals together, and 
this modification is necessary to ensure the mass distribu- 
tions of the HDE and LDE samples are comparable.) 

(vi) We assign all match galaxies to the LDE sample. If 
they were not already in the LDE sample, we set each of 
their weights to 1. Otherwise, we increase their weights by 
1. 

The resultant mass and redshift distributions of this 
scheme are assessed in Section [2.3.21 and Section [2.4.11 

2.3.2 Statistics of Galaxy Selection 

Fig. H] shows the distributions of stellar mass and redshift 
for the HDE and LDE samples of lens galaxies in the 
CFHTLenS. The matching scheme results in a nearly identi- 
cal distribution of stellar masses for HDE and LDE galaxies, 
and a very similar distribution of redshifts. 

Table[T]shows statistics for lens galaxies in the HDE and 
LDE samples in the CFHTLenS, for various stellar mass 
bins. The HDE sample contains a higher fraction of red 
galaxies than the LDE sample, as expected, but the dif- 
ference is at most 10% for a given stellar mass bin. This 
difference in the fractions of red and blue galaxies could in 
principle lead to a spurious detection of stripping if there is 
a relative bias in the stellar mass estimates between red and 
blue galaxies. This issue is discussed further in Section |4]4]T] 

2.3.3 Measuring the Lensing Signal 

To calculate the lensing signal around the HDE and LDE 
lens galaxies, we stack together all galaxies in a particular 
sample and stellar mass bir0. We then bin all lens-source 
pairs (only using pairs where Zphot.source > Zphot.iens + 0.1) 
based on the projected distance between the lens and source, 
calculated at the redshift of the lens. For each pair, we cal- 
culate the tangential ellipticity of the source relative to the 

* This process is performed one pointing at a time due to com- 
putational limitations, and all pointings are stacked together in 
the end. 
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SI.ellar Mass (M^„J Redshitt 

Figure 1. The distributions of stellar mass (left panel) and redshift (right panel) for the samples of HDE (solid line) and LDE (dashed 
line) galaxies, which, because of our matching algorithm, are virtually identical. The redshift distributions differ slightly between HDE 
and LDE galaxies, but there is no apparent trend to the deviation. 




Figure 2. The distributions of the types of galaxies classified as HDE and LDE in the simulations (left), and, of those classified as 
satellites, the distributions of the masses of the groups in which they reside (right). 



lens, gt, and convert this into units of surface mass density 
gradient: 

AE = Eeritf/t, (4) 
where 

E • - (5) 

and Da is the angular diameter distance to the source, Di 
is the angular diameter distance to the lens, and Dis is 
the angular diameter distance from the lens to the source. 



This me asurement relates to the projected mass of the lens 
through l|Mandelbaum et al.ll2005l ): 



(AE(i?)) = E(< R) - E(i?), (6) 



where E(< R) is the surface density averaged for all points 
contained within radius R, and E(i?) is the average surface 
density at radius R. This prescription works even for mass 
distributions that are not axisymmetric, as long as all points 
in a given annulus around a lens object are stacked. We 
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Table 2. Statistics of galaxies in the Millennium simulation for various stellar mass bins, using our models for estimating halo mass and 
environment, logm is the stellar mass bin. M is the mean halo mass of the galaxies in this bin in units of 10^*^ Mq, and z is their mean 
redshift. /sat, /fieldi £^nd /cen are the fractions of galaxies that are satellites, field galaxies, and group centrals, respectively. A/host is the 
mean mass of the host group for satellite galaxies in units of 10^^ Mq. All values assume galaxies are freighted by their stellar masses. 



compute the error in this value empirically from the scatter 
in calculated AE values for each lens-source pair. 

For the calculations of error in our model fits, we assume 
the noise in all radial annuli is independent. Strictly speak- 
ing, this isn't true, as there is a small correlation between 
the ellipticities of nearby sources, but this effect is negligible 
except at extremely large radial annuli. For computational 
simplicity, we do not apply the c2 correctiorjf] to source el- 
lipticities in our analysis. Because galaxy-galaxy measure- 
ments stack lens-source pairs over all position angles, they 
are insensitive to this correction(see IVelander et al.l |2012| . 
for further explanation and justification of this). Moreover, 
here we are interested in a differential measurement between 
galaxy-galaxy lensing samples, and so we expect our results 
to be highly robust to this effect. 



2.4 Simulated Galaxy Catalogues 

We require simulations in order to calibrate the frac- 
tions of satellite and central galaxies in our samples, and 
to test our methods for modelling the lensing signals 
around galaxies. T he simulations are b a sed o n the semi- 
analytic models of iDe Lucia fc Blaizoj (12007) which are 
based on the Millenium Simulation ( Springel et al.l l2005l : 
iLemson fc Virgo Consor tium 2 00^ ) . We use fo r our analysis 
a set o f thirty-two 16 deg'^ "lightcone" fields bv lHilbert et al.l 
(|2009l ). We assign photo-z errors consistent with those in 
CFHTLenS and apply the P3 algorithm to the simulated 
data. This allows us to select galaxies in the same manner 
as is done with the CFHTLenS dataset. 



2.4.1 Statistics of Simulated Catalogues 

Fig. [2] shows the distributions of galaxy types for the mock 
HDE and LDE samples drawn from the Millennium Simula- 
tion, and, for the satellite galaxies within each sample, the 
distribution of the masses of the groups in which they re- 
side. We classify galaxies as "central" (the most massive in 
a group), "satellite" (in a group but not the most massive) 
or "field" (not in a group). Table [5] shows the distributions 
of galaxy types for the HDE and LDE samples for various 
stellar mass bins. This shows that the fraction of satellites 



^ The c2 correction is an empirical correction to the e2 compo- 
nent of source ellipticity, based on the assumption that the mean 
e2 across a given field should be close to zero. 



in the HDE sample remains roughly constant with stellar 
mass, and decreases slightly with stellar mass in the LDE 
sample. For both samples, the fraction of centrals rises with 
stellar mass, while the fraction of field galaxies falls. Nei- 
ther sample shows any significant change with stellar mass 
in the mean mass of the host groups for satellites, except 
for a rise in the most massive stellar mass bin tested. HDE 
satellites are observed to reside in groups of ~ 4x 10" Mq. 
In contrast, for the small fraction of LDE galaxies that are 
satellites, the characteristic host halo mass is ~ 1.8 x 10^'^ 

Mq. 



2.4.2 Simulated Lensing Signals 

In order to simulate lensing signals for the Millen nium Simu- 
lation catalogues, we use the same methods as in lCillis et al.l 
120121). In short, we assume all galaxies and grou p centres are 



surrounded by spherical tr uncated NEW halos (IBartelmann 
1 19961: iHamana et al.ll2009l ). using the model of lBaltz et al 



( 2OO9I ) ■ and estimat ing halo m a ss fro m stellar mass by us- 
ing equation 3 from lCuo et al.l (|2010l ) (to better match the 
stellar-dark mass ratio in the CFHTLenS): 



0.129 X 



+ 



JTlhalo 



lOii-^M 



JTihalo 



1011-4M, 



(7) 



We form two versions of the simulated shear catalogue, 
one in which we simulate the effects of stripping by decreas- 
ing the truncation radii of satellite galaxies' halos and al- 
locating the lost mass to group centres' halos ("Stripping") 
and one in which we do not ( "No Stripping" ) . The amount of 
stripping is assumed to depend on distance fro m the centre 
of the host halo and is given by equation (5) of ICillis et al 



(2012), which is based on data in figure 15 of iGao et al 
(2004). With this prescription, we find that the mean re- 
tained mass after stripping is approximately 40% of the ini- 
tial mass. We then calculate shapes for all background galax- 
ies by assuming initially zero ellipticity in both components, 
then applying shear due to each nearby halo between the 
source galaxy and the observer. 
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Figure 3. An illustration of how the one-halo term varies with 
satellite mass (top left); and how the offset group halo term (equa- 
tion 1 II2I I) varies with group mass (top right), satellite concentra- 
tion (bottom left), and density threshold (bottom right). Plotted 
values of the parameters, with italicized parameter corresponding 
to the value used for other plots: Maat = 10^°, 10^^, 10^^ Mq; 

A/group = 10^'\ IQl*, 1015 Mq; c^at = 2.5, 5, 7.5, St = 0, 

10, 20 Mq/pc^. Increasing line weight corresponds to increasing 
the varying parameter. The fraction of satellites which reside in 
groups is not illustrated, as it is a simple scaling of the group halo 
term; it is fixed to 0.6 for these plots. 

3 MODELS AND SIGNAL FITTING 

We expect the lensing signal around galaxies in the HDE 
sample to be r easonably well-describ ed by the following halo 
model (see eg. IVelander et al]|2012l ): 

AE = AEih + /sat AEoG + AEah (8) 

where AEih is the "one-halo" term, /aat is the fraction of 
galaxies in the sample that are satellites AEqg is the "off- 
set group halo" term, and AE2h is the "two-halo" term, as 
described below: 

(i) One-halo term: The lensing signal that results from 
the galaxy's own dark matter halo. 

(ii) Offset group halo term: This is the contribution to 
the lensing signal a satellite caused by the presence of its 
group's halos. 

(iii) Two-halo term: Galaxies will typically reside near 
other massive structures, which results in a contribution to 
the lensing signal at large radii. 

Since galaxies in the HDE sample are more likely lie in 
overdense re gions, we cannot app ly exactly the same halo 
model as eg. iVelander et"ai] (|2012l l, who use all galaxies in- 
dependent of environment. This primarily affects the offset 
group halo term. See Section 13.21 below for an explanation 
of how we modify our halo model to account for this. 

For LDE galaxies, we expect the signal to be described 
by the form: 



where AEud is the "underdensity" term, which is the ef- 
fective contribution from the fact that galaxies in an under- 
dense environment will see a negative contribution to their 
lensing signal at large radii. This effect is analogous to the 
offset group halo term, except arising from an underdensity 
instead of an overdensity. 

We can best compare the lensing signals that result from 
stacks of HDE and LDE galaxies by fitting the signals with 
a model profile, and comparing these fits. The model profile 
for the HDE sample includes just the "one-halo" and "offset 
group halo" terms. Since the "underdensity" and "two-halo" 
terms are only significant at relatively large radii, we can 
safely ignore them if we do not fit the profiles out to large 
radii. We discuss this further in Section [4.4.21 

These components are discussed in the subsections be- 
low, and we discuss the procedure we use to fit a model to 
the data in Section [3.41 

3.1 One- halo term 

For the one-halo term, we assume that all galaxies reside in a 
dark matter halo that can be approxima ted with a tr u ncate d 
NEW density profile, as formulated by iBaltz et"al] l|2009l ). 
This model has three free parameters: the halo mass M200, 
concentration c, and the truncation parameter r = rtrunc/fg. 
In practice, we have found that the signal is not strong 
enough to simultaneously constrain all three parameters. 
Therefore, for simplicity, in our default fits discussed below, 
we fit M200, with c fixed by: 

taken from iNeto et al.l(|2007l ). and we also fix r = 2c, which 
is a reasonable value for unst ripped halos (|Hilbert fc Whit3 
I2OIOI : lOguri fc Hamanall201ll ) . 

In Section r4.3.2l we investigate alternative fits in which 
c or r are allowed to be depend on environment. 

3.2 Offset group halo term 

Since the P3 algorithm biases our galaxy selection such 
that the HDE sample predominantly consists of galaxies 
wit hin groups, we cann ot use the standard halo model 
(eg. IVelander et al]|2012l ) to calculate the contributions of 
nearby groups. Instead, we make the assumption that the 
sample consists of a fraction /aat satellites, and the rest are 
either central or field galaxies. The central and field galaxies 
will only have a one- halo component in their lensing signals, 
while satellites will have both the one-halo component and 
a contribution from their host groups. In order to model 
the average contribution of group halos to the lensing signal 
around galaxies in the HDE sample, we assume that it takes 
the following form: 

AEoG(i?) = AEi,,,tiR,Rs)P{Rs)dRs, (11) 

where Rs is the projected separation between a satellite and 
the group centre, AEhost(fi, ^s) is the contribution of the 
group halo to the lensing signal around a point at projected 
distance R from the group centre: 



AE = AEih + AEuD 



(9) 



AEhost (-R, Rs) — Ehoat (< R, Rs) — Ehost {R, Rs) 
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1 Z"^" 



2tt 



■^host 



(12) 



where Ehost(ii9) is the projected surface density of 
the host group's halo at projected radius Rg — 
VR'^ + Rs - R'Rs cos 9- and P{Rs ) is the probabihty that 
a sateUite in the sample will reside a distance Rs from its 
host group's centre. We assume P(-Rs) takes the form: 



Mgr, Csat)PHDE(-R3 



(13) 



where Mn is a normalization factor, E(i?s, Mgr, Csat) is the 
projected surface density of an NFW halo with mass equal 
to the mass of the host group, Afgr, but a concentration 
Caat, different from the dark matter concentration c. Analy- 
ses o f the satellite density in groups and clusters l|Lin et alj 
I2OO4I : iBudzvnski et al.1l2012l ) have indicated that the spatial 
distribution of satellites can be well-modelled in this way, 
assuming an NFW density profile with concentration ~ 2.5, 
which is lower than the typical concentration of the dark 
matter halo by a factor of ~2. 

The term Phdb(^?3) is the probability that a satellite at 
a distance Ra from the host group's centre will be included 
in the HDE sample. The form of PHDE(iis) is determined 
by the selection effects inherent in the P3 algorithm. To first 
order, P3 selects galaxies in regions of high projected surface 
density for the HDE sample. We thus model PHDE(i?s) as 
a smooth cut-off based on the projected surface density of 
the group. We wish for it to converge to PHDE(i?3) = 1 for 
T,{Rs) » T,t, and converge to PwEiRs) = for T,{Rs) ^ 
Et, and so we choose the following functional form, which 
has these properties: 



T.{Rs 



E(i?.)2 + E? 



(14) 



where T.{Rs) is the projected surface density for a satellite 
at distance Rs from a group centre and Et is the threshold 
density. As we have no prior justification for any specific 
density threshold to use, we leave this parameter free, to be 
fit by our algorithm. 

For the HDE sample, we fix /sat to the value found 
in the mock HDE sample from the Millennium Simulation. 
We do not expect this simulated result to perfectly match 
the fraction of satellites we might find in the CFHTLenS 
dataset, and we investigate what impact a different /sat 
might have in Section [4.4.21 For the LDE sample, we don't 
include this term, as the form of the measured lensing sig- 
nal in both simulated and CFHTLenS data shows that the 
underdensity signal dominates at large radii. 

We choose to model the offset group halo term as if 
all groups are of the same mass. We tested using a distri- 
bution of group masses, and the resulting signal was not 
appreciably different from the single-mass signal. The use of 
a distribution of group masses did tend to increase the re- 
sultant signal (the difference scaling with the spread of the 
mass distribution), even when the mean mass is fixed, and 
so the single-mass model will likely underestimate the mean 
host halo mass. 

Fig. [3] illustrates how the modeled one-halo term varies 
with satellite halo mass, as well as how the fitted offset group 



halo terms varies with the group halo mass Mgroup, satellite 
concentration Csat and threshold surface density Et. 

3.3 Underdensity signal 

Galaxies in the LDE sample are selected to lie in S/N < 
regions, which are underdense {S < 0) compared to a sur- 
rounding annulus with inner radius 1 Mpc and outer radius 
3 Mpc. Similarly to how galaxies in groups have a posi- 
tive contribution to their lensing signal from the offset over- 
density in which they reside, galaxies in underdense regions 
will have a negative contribution to their lensing signal on 
larger scales due to the fact that their local environment is 
less dense than the surrounding environment. This effect has 
been observed in both the CFHTLenS dataset, as well as in 
the simulations. 

The expected form of this negative lensing signal has 
not been well-studied, so there is no functional form which 
we expect it to take. We have attempted to fit this signal 
with the same functional form as the group halo term, mul- 
tiplied by a negative free term, but this failed to provide 
a suitable fit to either the simulated or to the CFHTLenS 
data. Note in the right panel of Fig. [5] that the minimum 
value for the LDE signal is at a higher projected radius than 
the peak of the offset group halo term. 

To handle this effect, for the LDE sample, we only fit 
the lensing signal for R < 400 kpc, where the one-halo term 
dominates the signal. 

3.4 Fitting Procedure 

For all fits, we use radial bins of 25 kpc< R < 2000 kpc. 
We tested constraining the fits to a lower maximum radius, 
and this had no noticeable effect on the fitted satellite halo 
masses. Fitting to a lower maximum radius only altered the 
fitted group mass, making it less well-constrained. 

We use a two-step procedure to fit the models to the 
lensing signals. Because our models are relatively simple, 
they are not perfect fits to the data. So, we first attempt 
to determine the amount of error inherent in our modeling, 
in order to assign more conservative uncertainties to the fit- 
ted parameters, as we will now describe. We first perform a 
steepest-descent minimization to obtain best-fitting pa- 
rameters for the model. At this point, if the Xrcd value for the 
fit is greater than 1, we assume that this is due to some error 
in the modeling, which we parametrize as am- We uniformly 
add this value in quadrature to the measured uncertainties 
in all radial bins, such that the adjusted x?ed = 1 for 
best fit. We then repeat this process, finding new best-fit 
values and recalculating cTm until convergence. 

Since this procedure effectively increases the error in all 
radial bins, this process has the result of increasing the mea- 
sured errors on all fitted parameters. If the model is initially 
a good fit (Xred ~ 1) to data, the increase is negligible, 
but if the model is a poor fit to the data, the estimated er- 
rors for the fitted parameters will be significantly increased. 
As such, this process allows us to place more conservative 
limits on our results, based on the quality of the model's fit 
to the data. 

Additionally, since the model error is uniformly added 
to the errors in all radial bins, it prevents the fitting algo- 
rithm from over-weighting the fit to the high-radius bins, 
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which otherwise have significantly lower errors, and thus 
typically contribute more to the x'^ value of the fit if the 
model isn't a perfect fit to the data. 

For the models we tested, we typically found for the 
HDE samples that am < 0.5Mq/pc^ , which is < 5% of the 
measured lensing signal AE. For the LDE samples, most fits 
were initially of Xred ~ li ^ind so no model error term was 
necessary. 

Once the model error is determined, we run an MCMC 
algorithm to help determine the errors of the fitted parame- 
ters. Since only the mass of satellite halos is relevant to us, 
we marginalize over all other parameters to get the mean 
value and errors for the satellite mass. 



4 RESULTS AND ANALYSIS 

In this section, we present the results of the fits and discuss 
their implications. In Section 14.11 we discuss the predicted 
results from the simulations, for both the "No Stripping" 
and "Stripping" models. Section 14.21 presents the main re- 
sults of our analysis of the CFHTLenS dataset and discuss 
their implications. In Section [4.31 we discuss alternative in- 
terpretations of the data, and which of the one-halo mass, 
concentration, and truncation radius might plausibly con- 
tribute to the observed differences between the HDE and 
LDE samples. Section [4.41 discusses potential systematic ef- 
fects. 



4.1 Predictions from Simulations 

Fig. [4] shows plots of the best-fit models for the simulated 
catalogues from the Millennium Simulation, for both the 
"Stripping" and "No Stripping models" (described in Sec- 
tion [2X2|) . for galaxies with W^Mq < m < 10^°-^Mq. The 
plot illustrates that in the "No Stripping" scenario, the mea- 
sured lensing signals for the HDE and LDE samples are 
nearly identical at very small radii. Our algorithm does not 
work perfectly for this mass bin, and in the "No Stripping" 
scenario, it fits a one-halo mass to the HDE sample that is 
somewhat larger than the one-halo mass fitted to the LDE 
sample, while for the "Stripping" scenario, the fitted HDE 
one-halo mass is slightly lower than the fitted LDE one-halo 
mass. 

Further comparisons of fitted one- halo masses for differ- 
ent mass bins can be seen in Table [3] and Fig. [T] As can be 
seen there, for all mass bins m < lO^^'^M© with the "Strip- 
ping" model, as expected the fit yields a relatively lower 
one-halo mass for the HDE sample compared to the LDE 
sample than it does for the "No Stripping" model. Above 
m = W^'^'^Mq, however, the fitted masses in the "Strip- 
ping" and "No Stripping" scenarios are comparable. This 
is due to the fact that at high stellar masses, the fraction 
of galaxies in the HDE sample that are centrals increases 
rapidly (see Table [2]). Since mass stripped from satellites is 
added to the masses of central galaxies, then if too many 
central galaxies are included in the sample, stripping will 
have little or no net effect on the lensing signal. 

The fitted group masses for the simulated data seen 
in Table [3] are larger than the actual group masses by a 
factor of ~ 1.5-2. Our tests have shown that this can occur 
when halos from a very broad range of masses are averaged 
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Figure 5. Measured lensing signal and model fits for data from 
the CHFTLenS survey, including all galaxies with 10^ < m < 



Mq. HDE (red) and LDE (blue) lensing signals and fits 



are illustrated, as well as the HDE data with the fitted offsct- 
group-halo term subtracted off (orange). The dashed line shows 
the one-halo model fit to the HDE sample, and the dotted line 
shows the HDE oflfset-group-halo term. The one-halo mass fit for 
the HDE sample is found to be significantly lower than for the 
LDE sample. 



together, as is the case here - the lensing signal of an average 
of halos of varying mass is similar to the lensing signal of a 
single halo with a mass somewhat greater than the average of 
the sample. The fitted group masses for the CFHTLenS data 
are additionally observed to be a factor of ~ 2 larger than 
the group masses for simulated data. This is not surprising, 
as the halo masses in the simulated data are extrapolated 
from the stellar masses of their constituent galaxies, and the 
distribution of stellar masses in the Millennium Simulation 
does not match the distribution in the CFHTLenS dataset. 

These results from the simulations imply that with the 
CFHTLenS data, a comparison of the HDE and LDE fitted 
one-halo masses can be used as an indication of whether or 
not tidal stripping is occurring, but we must use a stellar 



mass upper limit of ~ 10 M(? 



4.2 Observational Results 

Fig. [5] shows the lensing signals for the HDE and LDE sam- 
ples taken from the CFHTLenS survey, including all galax- 
ies with 10^ Mq <m < lO^"'^ Mq, with the best-fit models 
plotted on top. For this broad mass bin, the fits show that 
the HDE one-halo term is lower than the LDE term, at 2.5a 
significance (p = 0.0113). However, this simple fit is not op- 
timal. In part, this is because we are combining galaxies with 
greatly varying masses. The resultant lensing signal of this 
combination does not perfectly resemble the lensing signal 
of a single halo possessing the average mass of the sample, 
and the code compensates for this by fitting a higher a^i, 
which results in larger errors for the best fit. 
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Figure 4. Lensing signals and fits for simulated lensing data for the "No Stripping" (left) and "Stripping" (right) scenarios (see 
Section \2A\ . The "No Stripping" scenario shows similar one-halo fits for the HDE and LDE samples, while the "Stripping" scenario 
shows a lower one-halo mass fit for the HDE sample than for the LDE sample. Error bars are not shown, as shape noise is not simulated 
for these datasets, and so the scatter is extremely small. 



Table 3. Results of the fitting procedure when applied to simulated (top) and the CFHTLenS (bottom) lensing data in various stellar 
mass bins. All masses are in units of 10^" A^q. logm is the stellar mass bin. /sat is the fraction of satellites we use for the fitting, based 
on data from the Millennium Simulation. Mhde and Mlde are the fitted one-halo masses for the HDE and LDE samples. A/gr is the 
fitted mass of the offset group halo term. Rjvf is the ratio of AfjjDE to -'^^LDE- X^ed reduced parameter without the model error 

term (see Section l3.4p included (for 36 degrees of freedom; a value close to 1 is ideal). 
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Fig. |6] shows the Ukehhood distributions for the fit- 
ted satellite masses, host group mass, and surface density 
threshold for the HDE sample of galaxies with IQ^Mq < 
m < 10^"'^ M0. The plot shows that there is only a weak de- 
generacy of Maat with the other two parameters, but there is 
a stronger degeneracy between Mhoat and Et. Nevertheless, 
when marginalized over the other parameters, Mhost is very 
tightly constrained, and Msat is reasonably constrained. 

We can more carefully analyze the data by splitting 
the galaxy sample into smaller stellar mass bins. Fig. [7] 
shows the results of this analysis for both simulated and 



the CFHTLenS data, with the ratio of the fitted one-halo 
mass for the HDE sample to that of the LDE sample plotted 
against the galaxies' stellar masses. Simulated data is not 
available for all mass bins plotted due to limitations of the 
Millenium catalogue. The simulated data demonstrates that 
for sufficiently high stellar mass bins, the "Stripping" HDE 
mass becomes comparable to or greater than the "No Strip- 
ping" HDE mass. This is due to the fact that, at high stel- 
lar masses, the P3-identified HDE sample contains a large 
number of centrals. When tidal stripping is present, mass is 
transferred from satellites to centrals, increasing their mass. 



11 



10 
14.25 



A 






A 




(t5) 


^^^^ 


A 



1,4 11,5 11,6 11,7 10 12 14 16 14,1,5 14,20 14,25 

login Al,at.HDE (M.o,) El (Afc/pc^) logio A4o,t (Mg) 



Figure 6. Probability distribution functions and joint probability 
distribution functions for satellite mass Msat (left column), host 
group mass Mjjost (bottom row), and surface density threshold St 
(middle row and middle column) for the fit of the lensing signal 
of all HDE galaxies with lO^M© < m< 10^°-^Mq. 
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Figure 7. A summary plot of the fitting results for various stellar 
mass bins. Shown are the ratios of the best-fit one-halo mass 
for the HDE and LDE samples, for both the Stripping and No 
Stripping simulations, and for the CFHTLenS data. 



When centrals make up a large enough fraction of the sam- 
ple, the sample shows an increase in mass under the effects 
of tidal stripping. Therefore, we restrict the analysis to bins 
with m < lO^^'^Af©, where the "Stripping" scenario predicts 
lower fitted mass than the "No Stripping" scenario. 

Details of the fits to CFHTLenS data for different stellar 



mass bins are shown in Table [S] The goodness-of-fit is com- 
parable to previous galaxy-galaxy lensing studies. Specifi- 
cally, the Xicd values for our fits (which are calculated before 
the inclusion of the model error term, see Section 



f34ll 

similar to the full halo model fits of IVelander et all 



their x?ed values varied from 0.5-2 for different stellar mass 
bins, whereas ours vary from 0.6-2.3. 

If only the three stellar mass bins with 10^ < m < 
IO^^'^Mq are used, we obtained a weighted mean ratio of 
HDE one-halo mass to LDE one-halo mass of 0.65 ± 0.f2. 
If we assume that this ratio is indicative of the retained 
mass after stripping, and assume the sample contains ~ 60% 
satellites, then we can extrapolate that for a sample of 
100% satellites, the mean retained mass fraction will be 
~ 0.41 ± 0.19, which is consistent with the mean retained 
mass fraction of 0.40 we measured from the simulated data. 

Note that at face value, our result suggests less mass re- 
duction in HDE environments than the factors of 2-5 found 
for the ~ Lt satellites of the rich cluster CI 0024-1-16 found 
bv lNataraian et all Hooi). There are several key differences 
between these samples, however; in particular, our satellites 
have lower stellar mass and our satellites in habit lower mass 
host h aloes than the rich cluster studied bv lNataraian etld] 
(|2009l) . 

These combined results reject the results of the sim- 
ulated "No Stripping" model at 4.1(t (p < 0.0001), reject 
A^HDE = Mlde at 2.9(T (p = 0.0039), and are consistent 
with the simulated "Stripping" model at l.Scr {jp = 0.0651). 
This near-rejection of the "Stripping" model may indicate 
that this model underestimates the amount of tidal stripping 
which occurs in reality, or it might indicate that some effect 
other than tidal stripping (such as a difference in star for- 
mation histories dependant on environment) is contributing 
to the observed signal. Additionally, while our results have 
high statistical significance, they do not rule out the possi- 
bility of systematic errors resulting in a spurious detection. 
We investigate the possibility of such a spurious detection 
in Section [331 



4.3 Alternative Fits 

In the previous subsection, we made a number of assump- 
tions relating one-halo mass, concentration, and tidal radius 
(see equation H10|) ). This allowed us to compare only the one- 
halo mass between the HDE and LDE samples, to determine 
if there was a difference in their lensing signals. However, this 
doesn't tell us what variations in mass, concentration, and 
tidal radius might be causing this difference. In this section 
we investigate alternate means of fitting the lensing signals 
to determine which of these parameters might differ between 
the HDE and LDE samples. 



4-3.1 Truncation Radius 

An alternative fit with one free one-halo parameter is to fit 
the HDE model with the same one-halo mass as the LDE 
sample, but with a lower tidal (truncation) radius for its 
halo, as is predicted by typical models of tidal stripping. For 
the single broad stellar mass bin lO'^ < m < lO^''-^A/0 , we fit 
a tidal radius of rtidai/''20o = 0.26 ±0.14 for the HDE galax- 
ies. This corresponds to a typical tidal radius of 40 ± 21 kpc 
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Figure 8. Lower left: Joint probability distribution function for 
one-halo mass and concentration for the HDE (red solid lines) 
and LDE (blue dashed lines) galaxy samples. The dotted line 
shows the relation between mass and concentration given by equa- 
tion moi l, which was used for previous analysis. Upper left: Prob- 
ability distribution function for one-halo mass, marginalized over 
concentration. Lower right: Probability distribution function for 
concentration, marginalized over one-halo mass. 



and a retained mass fraction of 0.43 ± 0.18 for a satellite 
galaxy with M = 5.9 x 10^^ Mq and c = 8.5. However, this 
model is a marginally poorer fit compared to the default 
method of fitting one-halo mass: Xrcd = 1-08 for fitting tidal 
radius, compared to Xrcd = 0.99 for fitting mass. Our data is 
thus consistent with both interpretations, and we are unable 
to discern between them. 



4-3.2 One-halo Mass and Concentration 

It is also possible that the observed difference in lensing 
signals between HDE and LDE galaxies could be due to a 
difference in concentration, rather than a difference in mass. 
Halos with lower concentrations will have lower lensing sig- 
nals at small radii, and somewhat higher lensing signals at 
large radii. To test whether the observed results could be due 
to a change in concentration, we reran the analysis, leav- 
ing concentration as a free parameter. The resultant joint 
probability distribution function for this analysis is shown 
in Fig. [S] As this plot shows, while there is a difference in 
the mass PDFs for the HDE and LDE samples, there is 
no evidence for a difference in concentration between the 
two samples. We therefore favour the interpretation that the 
measured difference in lensing signals between the HDE and 
LDE samples is due to the HDE halos being less massive 
than LDE halos, but we are unable to rule out the possibil- 
ity that the two samples have the same mass, but the HDE 
sample is less concentrated. 



4.4 Analysis 

Although the results are statistically significant, it is nev- 
ertheless possible that systematic errors have entered our 
analysis. In this section, we discuss various possible system- 
atic errors that may affect the results. 



4-4-^ Stellar Mass Biases 

In principle, if there is a relative bias in the estimates of 
stellar mass between red and blue galaxies, and if the HDE 
and LDE samples contain different fractions of red and 
blue galax;ies, we could get a spuriously positive detection 
of tidal stripping with our method. It is difficult to com- 
pletely rule out the possibility of a relative bias in stellar 
mass, but it can be investigated by comparing the distribu- 
tions of stellar mass to other publis hed distributions. U nlike 
the CFHTLenS data, the WIRDS jSielbv et al.ll2012l ) data 
are based on optical photometry supplemented by deep in- 
frared images. We would therefore expect the stellar masses 
in WIRDS to be more accurate. A compariso n of stellar mass 
estimates between CFHTLenS and WIRDS (|Velander et al.l 
|2012| ) shows that, if we assume WIRDS stellar masses to be 
more accurate, then CFHTLenS stellar masses may indeed 
be slightly biased, with CFHTLenS red galaxy stellar masses 
~ 0.05-0.1 dex too low and the opposite for CFHTLenS 
blue galaxies. Since we are selecting by CFHTLenS stellar 
mass, and because red galaxies are more common in high- 
density environments, this implies that the true mean stellar 
mass in HDE regions is actually slightly larger than that in 
the "matched" LDE regions. Correspondingly, in absence of 
stripping, we would expect the recovered subhalo masses to 
be larger. Since we find them to be smaller, this effect is in 
the wrong sense to explain our stripping detection. 

In any case, the effect is small. As can be seen in Table[Tl 
the maximum difference in the fraction of red galaxies be- 
tween the HDE and LDE samples is in fact only 10%, in the 
lO''-'^ Mq < m < 10^° Mq and lO" Mq < m < 10^°-^ Mq 
bins. This, combined with the above estimates of the rela- 
tive bias, allows us to put an upper limit on this effect of 
-5%. 



4.4.2 Modeling Inaccuracies 

We have attempted to account for inaccuracies in the mod- 
els through the inclusion of a "model error" term, but this 
does not account for all possible errors in modeling that 
might arise. Notably, for the HDE sample, there is a weak 
degeneracy between the fitted one-halo term and the other 
fitted parameters, as well as with the fixed parameters: the 
fraction of galaxies which are satellites (/sat) and the con- 
centration of satellites (csat). We marginalize over the other 
fitted parameters to estimate the mean one-halo mass and 
its error, but errors in /sat and Csat would persist through 
the analysis as systematic errors. 

To assess the potential impact of errors in the fraction 
of satellites and satellite concentration, we reran the analysis 
varying these parameters, to observe how much the results 
changed. The results plotted here use a sample of all galax- 
ies with m < 1O^°-^M0, weighted by estimated halo mass. 
The results of this analysis are presented in Table |4l and 
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Figure 9. As Fig. [5] except with the fraction of sateUites fixed to 0.4 (left), and the concentration of satellites fixed to 5.0 (right) for 
the fitting of the HDE signal. 



Table 4. Results of varying the fixed parameters (/sat and Csat) 
in the fitting procedure, using a sample of galaxies with 10^ A/0 < 
m < lO^"'^ Mq. All masses are in units of lO"' Mq. 



fsat 


Csat 






A/group 




0.4 


2.5 


44 ± 6 


0.76 


23000 


1.17 


0.6 


2.5 


38 ± 5 


0.65 


15000 


0.99 


0.8 


2.5 


33 ± 5 


0.56 


12000 


1.22 


0.4 


5.0 


35 ± 4 


0.59 


25000 


1.08 


0.6 


5.0 


25 ± 4 


0.43 


17000 


1.93 


0.8 


5.0 


20 ± 4 


0.33 


9100 


12.09 



Table 5. Fitted values of St for various stellar mass bins, in units 



of Afg/pc^. 




log m 


St 


9-9.5 


42 ± 2 


9.5-10 


34 ± 3 


10-10.5 


44 ± 3 


10.5-11 


39 ± 6 


11-11.5 


36 ± 7 



sample plots of the fits can be observed in Fig. |9l We al- 
low tlie fraction of satellites to vary from 0.4 to 0.8, and we 
test increasing the concentration of satellites to 5.0, which 
is comparable to the concentration of the dark matter halos 
of galaxy groups. 

As compared with the default case (/sat = 0.6 and 
Csat ~ 2.5), the only variation of the parameters that re- 
sults in increasing the fitted HDE mass is when the fraction 
of satellites is decreased. If the tested fraction of 0.4 were 
the case in reality, the significance of the detection would 
decrease to ~ 2.1a. However, this particular satellite frac- 
tion results in a marginally poorer fit, as evidenced by the 
increased x?cd value. 

Alternatively, if the actual concentration of satellites 
within groups is higher, the fitted satellite mass will be 
lower. A fit with Csat = 5.0 yields a marginally poorer Xrcd 
value, but we cannot rule out this scenario. If this is the 
case in reality, it would strengthen the significance of the 
detection to ~ 4.4(t. 

It is also important to look at the impact of the den- 
sity threshold term, Et. Although this term was found to be 
necessary in the simulated data to provide a reasonable fit 
when the fraction of satellites is known exactly, it is not cer- 
tain that this is an accurate description. We can investigate 



this matter by looking at how the fitted value of Et varies in 
the fits, as can be seen in Table [5] If we were modeling ev- 
erything perfectly, we would expect to see Et being roughly 
constant for all stellar mass bins, as the P3 algorithm is blind 
to galaxy mass, and so there is no reason the threshold Et 
for S/N > 2 regions should vary with the mass of satel- 
lites. In fact, we observe Et to remain roughly constant at 
~ 4OM0/pc^, which is consistent with the hypothesis that 
we are modeling it reasonably. 

Other modeling inaccuracies may also affect the results. 
The model for the HDE sample neglects the contribution of 
the two-halo term (the contribution of nearby groups and 
field galaxies) to the lensing signal around satellites. We 
tested the implications of this with a rough model of the 
two-halo term, and it resulted in decreasing the fitted HDE 
satellite mass. Therefore, any possible systematic error in 
our results from this effect would be in the wrong sense to 
contribute to a spurious detection. 

Additionally, the model for the LDE sample neglects the 
contribution of the local underdensity (see Section I3.3|l to 
the lensing signal. Proper handling of this term would likely 
result in a slight increase in the fitted LDE mass, which 
would increase the significance of the detection. 
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5 CONCLUSIONS 

Previous lensing analyses of t he environmental depen- 
dence of satellite halo masses l|Mandelbaum et alj l2006l : 
Ivan Uitert et al.1l201lh have revealed the difficulty of detect- 
ing such an effect in samples where the satellites are predom- 
inantly expected to lie in groups. The analysis here improves 
on previous work by using the much deeper data provided 
by the CFHTLenS sample to better constrain the lensing 
signal around low-mass satellite galaxies in high-density en- 
vironments. 

Using photometric redshifts we divide galaxies in high- 
density and low-density environment subsamples that are 
matched in stellar mass. We have found a significant dif- 
ference in their halo masses. Our analysis shows a highly 
significant (AAa) rejection of the simulated "No Stripping" 
model, and a significant (3.2(t) rejection of the simple "null 
hypothesis" that there is no difference in the halo proper- 
ties of IfDE and LDE galaxies, for galaxies in a broad range 
of group masses. This difference is most likely due to tidal 
stripping of dark matter, and if so, this analysis represents 
the first detection of tidal stripping in a selection of galaxies 
that do not all reside within galaxy clusters. 

We argue that these results are unlikely to be due to 
systematic errors in our methodology, as most suspected sys- 
tematic errors would tend to bias one against a detection of 
stripping. 

The mean ratio of fitted mass for the high-density en- 
vironment sample to that of the low-density environment 
sample is ~ 0.65 ± 0.12. Since the HDE galaxy sample con- 
sists of only ~ 60% satellites, the retained mass fraction 
for a pure satellite sample would be considerably lower: 
~ 0.41 ± 0.19. We can alternatively model this as the 
HDE satellites being tidally stripped at a typical radius 
of rtidai/r200 = 0.26 ± 0.14. This corresponds to a typical 
tidal radius of 40 ± 21 kpc and a retained mass fraction of 
0.43 ± 0.18 for a satellite galaxy with M = 5.9 x 10"Mq 
and c = 8.5. 

Further work will be necessary to confirm these results, 
and to analyze the dependence of stripping on both group 
mass and on the satellite's location within the group. 
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